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ABSTRACT 

We report the detection of a Jupiter-mass planet discovered with the SOPHIE spectrograph mounted on the 1.93-m telescope at the 
Haute-Provence Observatory. The new planet orbits HD 109 246, a GOV star slightly more metallic than the Sun. HD 109 246b has 
a minimum mass of 0.77 Mjup, an orbital period of 68 days, and an eccentricity of 0.12. It is placed in a sparsely populated region 
of the period distribution of extrasolar planets. We also present a correction method for the so-called seeing effect that affects the 
SOPHIE radial velocities. We complement this discovery announcement with a description of some calibrations that are implemented 
in the SOPHIE automatic reduction pipeline. These calibrations allow the derivation of the photon-noise radial velocity uncertainty 
and some useful stellar properties (v sin /, [Fe/H], log -KJj,^) directly from the SOPHIE data. 
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1. Introduction 

I Fifteen years after the first discovery of a planet orbiting another 
solar-type star (Mayor & Queloz 1995), about 450 extrasolar 
planets have been announced. Most of them were detected with 

, the same historical technique: high-precision radial velocime- 
try. Although several methods now allow the routine detection 
of extrasolar planets, high-precision radial velocimetry will re- 

, main at the forefront of exoplanet science in the coming years. 
Indeed, due to continuous efforts to improve the instrumental 
sensitivity (e.g. Mayor et al. 2009) and to overcome the limita- 
tion of stellar activity (e.g. Udry et al. 2006), radial velocimetry 
is a particularly efficient technique to detect low-mass planets 

' (e.g. Mayor et al. 2009). Furthermore, radial-velocity (RV) mea- 
surements are needed to establish the true nature of transiting 
planetary candidates, to derive their mass, and to measure spin- 
orbit alignment angles via the Rossiter-MacLaughlin effect (e.g. 
Queloz et al. 2000; Ohta et al. 2005; Winn 2010; Triaud et al. 
2010). The numerous projects of future high-resolution spectro- 
graphs operating either at infrared (SPIRou, UPF, SIMPLE) or 
at visible wavelengths (ESPRESSO, CODEX) denote the field 
development and illustrate the importance of RV measurements 
in exoplanetary science. 

* Based on observations made with the SOPHIE spectrograph on the 
1.93-m telescope at Observatoire de Haute-Provence (CNRS/OAMP), 
France (program 07A.PNPCONS). 



The high-precision SOPHIE spectrograph (Perruchot et al. 
2008; Bouchy et al. 2009b) has replaced ELODIE (Baranne et al. 
1996; Queloz et al. 1998b) at the Cassegrain focus of the 1.93- 
m telescope at the Haute-Provence Observatory (OHP, France). 
Opened to the community since October 2006, SOPHIE has 
led to many discoveries in various research fields, including the 
follow-up of transiting planet candidates from the SuperWASP 
(e.g. Cameron et al. 2007), HAT (Bakos et al. 2008) and CoRoT 
(e.g. Bouchy et al. 2008) surveys; asteroseismology (Mosser et 
al. 2008); and the dynamics of binary stars (Albrecht et al. 2009). 
In late-2006, the SOPHIE consortium initiated a large program 
to search for and characterize extrasolar planets (Bouchy et al. 
2009b). This comprehensive survey is both a continuation and 
an extension of the planet-search programs carried out with the 
ELODIE spectrograph (e.g. Mayor & Queloz 1995; Delfosse et 
al. 1998; Perrier et al. 2003; Galland et al. 2005; Naef et al. 2005; 
Da Silva et al. 2006). 

HD 109246 was observed as part of the second subprogram 
from the SOPHIE consortium. This subprogram comprises two 
parts: a survey for giant planets in a volume-limited sample of 
~2000 FGK dwarfs, and a follow-up of known transiting gi- 
ant planets (Bouchy et al. 2009b). Previous discoveries from 
the main survey include the detection of a massive planet (or 
light brown dwarf) around HD 16760 (Bouchy et al. 2009b) and 
the detection of a two-planet system around HD 9446 (Hebrard 
et al. 2010a). The follow-up part allowed the observation of 
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the spectroscopic transits (Rossiter-McLaughlin effect) of three 
planets and revealed the two first cases of spin-orbit misaligne- 
ment (Loeillet et al. 2008; Hebrard et al. 2008; Moutou et al. 
2009; Pont et al. 2009; Hebrard et al. 2010b). 

Here, we report the discovery of a Jupiter-mass planet orbit- 
ing HD 109246. The observations are presented in Sect. 2 and 
the characteristics of the host star are described in Sect. 3. In 
Sect. 4 we derive the best-fit Keplerian model to the velocities of 
HD 109246 and correct for an instrumental effect on radial ve- 
locities due to seeing variation. We summarize and discuss our 
results in Sect. 5. In the Appendices, we describe the calibra- 
tions implemented in the SOPHIE automatic reduction pipeline 
to derive the photon-noise RV uncertainty ([A]), the projected ro- 
tational velocity v sin / ( IB. lb . the metallicity index [Fe/H] ( IB.2l i. 
and the activity index logR'^^ Q. 

2. Radial-velocity measurements 

We observed HD 109246 with the SOPHIE spectrograph from 
January 2007 to February 2010. The measurements were per- 
formed in high-resolution mode (resolution power of AA/A » 
75 000) and at constant signal-to-noise ratio (SNR variations 
CsA'R-8%). Exposure times ranged from 270 to 1200s, follow- 
ing the variations of seeing and sky transparency at OHP We 
used the objAB observing mode with one of SOPHIE'S fiber 
recording the starlight and the other fiber observing the sky to 
estimate the background contamination (e.g. moonlight). Both 
fibers were fed with a Thorium- Argon lamp every ~ 2 hour for 
wavelength calibration. SOPHIE is environmentally stabilized 
and the intrinsic drift of the spectrograph (mainly due to small 
variations of pressure and temperature) is less than 3 ms"' per 
hour. Each exposure was corrected for the instrumental drift by 
interpolating the measured intrinsic drift over the ~ 2 hour pe- 
riod. This method allows a precision of ~2ms ' on the correc- 
tion of the instrumental drift as determined empirically compar- 
ing the interpolated drift with the measured one from observa- 
tions made with simultaneous wavelength calibration. 

We used the SOPHIE automatic pipeline (Bouchy et al. 
2009b) to extract the spectra from the detector images and to 
cross-correlate them with a G2-type mask derived from the Sun 
spectra. We obtained the radial velocities by fitting each result- 
ing cross-correlation function (CCF) with a Gaussian (Baranne 
et al. 1996; Pepe et al. 2002). At this stage we rejected three 
spectra. Two of these spectra had a low signal-to-noise ratio 
(SNR(550nm)<30), which rendered the correction for the CCD 
charge transfer inefficiency (CTI) less effective (Bouchy et al. 
2009a). The third spectrum was contaminated by moonlight. 

The final RV data set for HD 109 246 comprises 58 measure- 
ments with a typical SNR of 47 (per pixel at 550 nm). This leads 
to a mean measurement uncertainty of 3.9ms"', including pho- 
ton noise and wavelength calibration errors. The photon-noise 
uncertainty was estimated as described in AppendixlAl External 
systematic errors of 2 ms"' (spectrograph drift uncertainty) and 
of 4 ms"' (guiding errors) were quadratically added to the mean 
measurement uncertainty (Boisse et al. 2010a). The mean error 
bar on each RV measurement is then about 6 ms"' . The RV mea- 
surements of HD 109 246 are listed in Table IC. 11 available at the 
CDS. Table lC. ll contains in its cols. 1-3, the time of the observa- 
tion (barycentric Julian date), the RV and its error, respectively. 

3. Stellar properties 

HD 109 246 (HIP 61 177) is a GO dwarf with an apparent V-band 
magnitude of my = 8.77 (SIMBAD database) and an astromet- 



Table 1. Stellar parameters for HD 109246. (I): Parameter de- 
rived from the SOPHIE CCF (see Appendices). (2): Parameter 
derived from the SOPHIE spectra (see Appendices). 



Parameters 


Values 


Sp.T 


GOV 


mv 


8.77 


B - V 


0.64 


7T [mas] 


15.24 ±0.68 


T.// [K] 


5844 ±21 (2) 


log g [cgs] 


4.46±0.19<2' 


Fe/H [dex] 


0.10 ±0.05 


V sini [kms"'] 


3±1(') 


M. [Mo] 


1.01±0.11 


R. [Ro] 


1.02 ±0.07 


log^HK [dex] 


-5 .05 ±0.1 


Distance [pc] 


65.6 



ric parallax of ;r = 15. 14 ± 0.68 mas measured by Hipparcos (van 
Leeuwen 2007). Allende Prieto & Lambert (1999) derived a stel- 
lar mass ]VI. = 1 .01 ± 0. 1 1 IVIq and a radius R. = 1 .02 + 0.07 Rq. 

We determined the star's effective temperature, grav- 
ity and metallicity using the spectroscopic analysis of 
Santos et al. (2004). This analysis was performed on the 
high-SNR average spectrum obtained by summing all the 
SOPHIE spectra. The spectroscopic analysis gave the fol- 
lowing results: an effective temperature Tgff = 5844 ± 21 K, 
a surface gravity log g =4.46 + 0.19, a micro-turbulence 
velocity V, = 1.01 ±0.02 km s"', and a stellar metallicity 
[Fe/H] =0.10 + 0.05 dex. When combined with isochrones (da 
Silva et al. 2006fl these parameters yield a stellar mass 
M, = 1.04 + 0.10 Mq, in agreement with the Allende Prieto & 
Lambert (1999) value. We also computed the projected rota- 
tional velocity and the stellar metallicity from the SOPHIE CCF 
as described in Appendix iB] This gave v sin / = 3 ± 1 km s"' and 
[Fe/H] = 0.14 ± 0.09. This alternative estimation of the stellar 
metallicity is consistent with the accurate determination based 
on spectral analysis. 

We estimated the stellar activity level from the emission in 
the core of the Call H&K bands measured on each spectra 
of HD 109246. The calibration used to derive logR'^^^ activ- 
ity indexes from SOPHIE spectra is described in Appendix ICl 
Applied to HD 109246, this calibration yields log/?[jj,=- 
5.05 + 0. 1, which translates into an RV dispersion that would be 
smaller than 1 ms"' according to our experience on HARPS. The 
stellar parameters adopted for HD 109 246 are listed in Table [T] 

4. HD 109246b's orbital solution 

4.1. Evidence for a planetary companion 

As shown in Fig. 1, the RV data of HD 109 246 exhibit a peak-to- 
peak variation of ~95 m s"' and a dispersion of ctrv =28 m s"' . 
A clear signal with a semi-amplitude larger than 30 ms"' is 
identified in the RV periodogram around 68 days. The disper- 
sion of the bisector span values (BIS) is relatively low, with 
(Tbis =9.9 ms"' (Fig.|2]i. No periodicity is detected in the bisec- 
tor periodogram and no correlation is present in the bisector-RV 
plot (Fig. |2]l. All these observations support a scenario where 
the RV variations of HD 109 246 are caused by the gravitational 
perturbation of an orbiting companion. 

We fitted the RV data of HD 109 246 with a Keplerian model 
using a Levenberg-Marquardt algorithm, after selected starting 

' Web interface available on ' http://stev.oapd.inaf.it/cgi-bin/paraml 
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Fig. 1. Radial velocities (top) and residuals to the best-fit 
Keplerian model {bottom) for HD 109246 as a function of 
barycentric Julian date. The best-fit Keplerian model is repre- 
sented by the black curve. The residuals have a dispersion of 
10.2ms-i. 

Table 2. Keplerian solution and inferred planetary parameters. 
Results in the first column are based on the original RVs. Results 
in the second column are based on the corrected RVs (see text 
for details). 



Parameters 


Original velocities 


Corrected velocities 


VR„ [kms"'] 


-19.463 ±0.004 


-19.464 ± 0.002 


P [days] 


68.20 ± 0.17 


68.27 ±0.13 


K [ms-'] 


38.2 ± 2.2 


38.2 ± 1.6 


e 


0.17 + 0.07 


0.12 ±0.04 


oj [deg] 


273 ± 25 


235 ± 29 


To [JD] 


54831.7 ±4.2 


54824.6 ± 4.7 


Mp sin / 1 [Mj„p] 


0.76 ±0.11 


0.77 ± 0.09 


a' [AU] 


0.33 ± 0.09 


0.33 ± 0.08 


a-(o-o [ms"'] 


10.2 


7.7 



assuming M. = 1 .0 1 ± 0. 1 1 Me 



values with a genetic algorithm. The best solution is a mildly 
eccentric orbit (e = 0.17 + 0.07) with a period P = 68.20 + 0.17 
days and a semi-amplitude K = 38.2 ± 2.2 ms"'. This signal cor- 
responds to a planet of minimum mass Mp = 0.76 ± 0. 1 1 Mjup or- 
biting HD 109 246 with a semi-major axis of 0.33+ 0.09 AU. The 
best-fit Keplerian model is plotted superimposed to the SOPHIE 
velocities in Fig.[T] The corresponding phase-folded RV curve is 
shown in Fig.|3]and the orbital elements of the planet are listed in 
Table |2] Uncertainties correspond to the 0.95 confidence interval 
after 5000 Monte-Carlo simulations. The tool used for the model 
fitting was successfully employed in exoplanetary surveys (e.g. 
Mayor et al. 2009; Bouchy et al. 2009b; Hebrard et al. 2010a). 
The uncertainties in the minimum mass and in the semi-major 
axis take into account the uncertainties both in the stellar mass 
and in the signal semi-amplitude, period and eccentricity. 

The residuals to the best-fit Keplerian model, o-(0 - C) - 
10.2ms"', are a bit large compared to the mean error bar of 
«6ms"'. The Tednced is equal to 1.8, which suggests that 
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Fig. 2. Bisector span (BIS) as a function of RV (top), and as a 
function of the RV residuals (bottom). Uncertainties in the BIS 
are assumed to be twice the RV ones. In both panels the scale is 
the same in the x and y axes. No correlation is seen between BIS 
and RV residuals. The extra uncertainty seen in the RV data is 
unlikely to be caused by stellar activity. 



an additional variation with a dispersion of ~8ms"' is present 
in the data. If stellar jitter was the source of the excess vari- 
ability (in spite of an estimated value lower than 1ms"', see 
Sect. [3]), we could expect a negative correlation between the bi- 
sector span and the RV residuals (e.g. Queloz et al. 2001; Boisse 
et al. 2010b). As shown in Fig.|2] no such correlation is observed 
in the data. We note that the expected correlation would be very 
close to the detection limits but the level of activity needed to 
cause the RV variations is inconsistent with the logR'^^ value. 

About a quarter of gaseous giant planet are in a multi- 
ple planetary system. The excess variability may thus alterna- 
tively be due to the presence of an additional planet in the sys- 
tem. We did not find any indication for a second planet in the 
present data set. With a maximum semi-amplitude of 14ms"', 
the RV residuals exclude the presence of an inner planet with 
mp sin i > 0.28 Mjup. A longer period planet could not induce a 
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Fig. 3. Phase-folded RV curve for HD 109246. The soHd Hne 
represents the Keplerian fit with a reduced equal to 1.8. The 
period is 68.2 d, the eccentricity is 0.17 and the semi-amplitude 
is 38.2 m s"' . Individual error bars are also plotted with a typical 
value of 6ms 
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Atmospheric condition estimator S 

Fig. 4. Residuals to the Keplerian fit as a function of the atmo- 
spheric condition estimator S, inversely proportional to seeing 
(see Sect. l4.2l for the definition of 2). The least-square fit repre- 
sented by the solid line has a correlation coefficient of -0.6 and 
a Spearman coeficient of -0.56. 

linear drift larger than 0.4 ms Additional data will be needed 
to probe the presence of low-mass inner planets and to constrain 
the presence of possible outer companions. 

The third and best possible explanation to the relatively large 
RV residuals is the underestimation of an instrumental effect. 
Described in the next Section, the so-called seeing effect seems 
to affect the observations of HD 109 246 (Fig. |4]l. 

4.2. Spectrograph illumination limitation 

The large residuals of HD 109 246 RVs to the best-fit Keplerian 
model seem well explained by seeing effect due to bad scram- 
bling of one multimode fiber (Fig.|4|. This effect was presented 
in details in Boisse et al. (2010a). We summarize here the main 
issue for the SOPHIE HR mode. 

SOPHIE is a fiber-fed spectrograph. The stellar light col- 
lected by the telescope is led to the instrument through a standard 
step-index multi-mode cylindrical optical fiber. They have good 



azimuthal scrambling but the radial one is not perfect. In order 
to improve the scrambling ability of the fiber and to stabilize the 
illumination, optical double scrambler are usually coupled to the 
fibers (as it is done for the SOPHIE HR mode) (Brown 1990; 
Hunter & Ramsey 1992). Non-uniform illumination of the slit 
or output fiber at the spectrograph entrance (variations in seeing, 
focus and image shape) decreases the radial-velocity precision, 
due to non-uniform illumination inside the pupil of the spec- 
trograph that lead, with optical aberrations to variations in the 
centroids of the stellar lines on the focal plane. 

Our experience on SOPHIE has led to identify remaining RV 
limitations due to the incomplete fiber scrambling. We simulated 
the optical path in the SOPHIE spectrograph and observed that 
variations of the slit, pupil or optical fiber illumination are di- 
rectly translated on the spectrum. The displacement of the spec- 
trum on the CCD detector is function of the wavelength (or the 
position on the detector) and is proportional to the width of the 
input image at the entrance of the fiber. Because the effect is not 
symmetric along the spectral orders and not monitored by the 
calibration lamp, the final computed RV varies and its variation 
is proportional to the input image width. 

The value of the seeing is not monitored by SOPHIE up until 
recently. We estimated the atmospheric condition by calculating 
the relative flux by unit of exposure time: 



r„pio-«''/2.5 



, or for the same star Z = 



SNR^ 
T 



(1) 



with SNR, the signal to noise ratio of the spectra, T^xp the ex- 
posure time of the measurement. My the visual magnitude of 
the target. This parameter allows a relative estimation of the see- 
ing since we cannot evaluate the atmospheric absorption. The 
atmospheric condition estimator 2 is inversely proportional to 
the seeing value. 

The seeing effect is expected to affect the measurements 
since the input beam at the entrance of SOPHIE is center- 
illuminated (input image smaller than the 3" sky acceptance of 
the optical fiber). In this case, simulations predict a systematic 
decrease in measured RV values, which is observed in the ob- 
tained data, as shown in Fig.|4]of this paper and for a sample of 
SOPHIE targets in Fig. 5 from Boisse et al. (2010a). Moreover, 
when the input image is smaller than the diameter of the fiber, we 
are more sensitive to effect of guiding and centering system that 
introduce RV variations. These RV variations are quite random 
and add some noise, whereas the seeing effect is directly related 
to the optical path in the spectrograph and may be quantified. We 
are modeling these variations with our data in order to remove 
this noise with a software tool even if this correction is only rela- 
tive due to our inability to estimate the absorption parameter. The 
seeing effect may reach peak-to-peak RV variations of 20 ms"' 
and it is the main current limitation of the SOPHIE accuracy. 
A new guiding camera is now in operation on SOPHIE and al- 
lows a better guiding and centering and also a monitoring of the 
exact seeing. It will then allow us to refine this method of cor- 
rection. Tests on square and octagonal section fibers to optimize 
the fibers scrambling are now under development and SOPHIE 
will be used as a bench test to validate these new feed optics. 
Soon, a new double scrambler will be mounted on SOPHIE with 
the goal to reduce this effect by at least a factor 10. 

4.3. Corrected velocities and improved orbital solution 

In Fig. 2] the residuals to the Keplerian fit are correlated with 
the atmospheric condition estimator 2 (see Eq. [1]). We calcu- 
late the linear correlation coefficient that is equal to -0.6, and 
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Fig. 5. Mass limit exclusion diagram with 99% confidence level 
for an additional companion in the HD 109 246 system calcu- 
lated from the residuals after the best-fit Keplerian model on the 
corrected velocities. The dashed line corresponds to the obser- 
vation span. The two days peak expresses the lower efficiency at 
this period due to data sampling. 
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Fig. 6. Period distribution of the known extrasolar planet orbit- 
ing single dwarf stars discovered via RV surveys. With a period 
of 68 days, HD 109246b is in the middle of the "period valley". 
Horizontal stripped illustrates the distribution for planets more 
massive than 0.8 M y„p and the shaded one for planets more mas- 
sive than 4 Mjup. 



the Spearman coefficient equal to -0.56, that support the anti- 
correlation between the two parameters. The linear function de- 
termined by a least-squares fit was subtracted from the radial- 
velocity data measurements to correct for the seeing effect. 
Fitting a Keplerian model to the corrected RVs yields the orbital 
parameters reported in the right column of Table 3. The new 
planetary orbit is similar to the previous one; the major change 
introduced by the correction is to decrease the RV residuals to 
7.7 ms"' . This new value does not change a lot our previous con- 
clusions regarding the possible presence of an additional planet 
as illustrated in Fig. |5] The remaining variability excess might 
be accounted for by an increase in guiding and centering er- 
rors when the input image is smaller than the fiber diameter (see 
Sect. 14.2b or by RV jitter due to stellar activity (see Sect.[3]l. 

5. Discussion and conclusion 

We report the detection of an extrasolar planet candidate around 
HD 109 246 discovered with the SOPHIE spectrograph mounted 
on the 1.93-m telescope at the Haute-Provence Observatory. 
HD 109 246b is characterized by a minimum mass of 0.77 ± 0.09 
Mjup, an orbital period of 68.27 ±0.13 days, and an eccentric- 
ity of 0.12 + 0.04. The host star being slightly more metallic 
than the Sun, the discovery of HD 109 246b reinforces the cor- 
relation between giant planet occurrence and stellar metallicity 
(e.g., Santos et al. 2004; Santos et al. 2005; Fischer & Valenti 
2005). No photometric search for transits have been undertaken 
until now though the transit probability is slightly above 1%. 

The residuals to the single-planet Keplerian fit (10.2ms ') 
are larger than the mean error bar of « 6 m s"', which suggests 
the presence of an additional source of RV variability. We con- 
sidered three possibilities: stellar jitter, the presence of an addi- 
tional planet in the system, and an instrumental effect related to 
seeing variations. We did not detect a clear signature of stellar 
activity, nor that of an additional planet. On the other hand, we 
identified a correlation between the residuals to the Keplerian fit 
and a seeing indicator. Correcting the RVs of HD 109 246 for this 
systematic effect decreased the RV residuals to 7.7ms"' while 



preserving the original planetary orbit. This result indicates that 
the seeing effect described in Boisse et al. (2010a) is at least 
partly responsible for the excess variability measured in the ve- 
locities of HD 109 246. A new guiding camera is now in oper- 
ation on the 1.93-m telescope at OHP Besides providing bet- 
ter guiding and centering performances, this new camera allows 
for a direct monitoring of the seeing. Future observations shall 
thus allow the method described in this paper to correct SOPHIE 
data for the seeing effect to be refined. Moreover, a new double 
scrambler will be mounted on SOPHIE to improve the scram- 
bling ability of the fibers and remove the seeing effect. 

The most remarkable property of HD 109 246b is that its or- 
bital parameters place it within the so-called "period valley" in 
the period distribution of extrasolar planets (Fig.|6l). The period 
valley extends from ~10 to ~100 days (Udry et al. 2003). It sep- 
arates the pile-up of hot Jupiters from the population of giant 
planets with periods above ~1 year HD 109246b mass rein- 
forces the fact that light planets (mp sini < 0.8 Mjup) are pref- 
erentially found to orbit with short period (P < lOOd) (Udry 
& Santos 2007). Ignoring multiple stars systems, RV surveys do 
not detect massive planets with nip sin i > 4 Mjup with period 
less than 100 days. The two components left are HD 162020b, 
a probable brown dwarf (Udry et al. 2002) and one of the two 
massive components orbiting HD 168443 (Marcy et al. 2001; 
Udry et al. 2002). Including the transit surveys detections, plan- 
ets more massive than 4Mjup are not found in the "period val- 
ley". 

A comparison with the models of Mordasini et al. (2009) 
shows that HD 109 246b may belong to the population of "main 
clump" planets. If so, it probably started to form near 4-6 AU 
and then moved inwards through type II migration. As outlined 
by Mordasini et al. (2009), the number of giant planets orbit- 
ing inside 1 AU is an important observational parameter for 
planet formation models. Additional detections of planets like 
HD 109246b will thus help characterize the profile of the sur- 
face density of solids below 1 AU, and the efficiency with which 
migrating cores can accrete planetesimals. 
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Appendix A: Determination of tlie RV uncertainty 
from the SOPHIE Cross-Correlation Function 

The methodology to compute the photon-noise uncertainty of 
radial-velocity measurement was described by Bouchy et al. 
(2001). In this approach, the quality factor Q represents the in- 
strinsic RV-information content of a given spectrum (quality and 
richness of spectral lines). This Q factor is independent of the 
flux and is computed on a reference spectrum considered as 
noise free. The RV uncertainty 6Vrms is then computed using 
the following relation: 



6V, 



RMS 



(A.l) 



with c the speed of light and Ne- the total number of photoelec- 
trons counted over the whole spectral range. In practice one does 
not always have a noise free reference spectra and one needs to 
estimate correctly the RV uncertainty on real observed spectra. 
On a noisy spectrum, the high frequency structures due to noise 
increase artificially the Q factor and lead to an underestimated 
uncertainty for the RV. We estimate that the computation of the 
Q factor is correct for spectra with S/N ratio per pixel higher than 
~200. For lower SNR, Q is overestimated and SVrms is therefore 
underestimated. 

We adapted our methodology directly to the Cross- 
Correlation Function which corresponds to an average spectral 
line in velocity space at much higher S/N than individual lines 
and may be considered as noise free for the computation of the 
quality factor. In that case Qccf can be expressed by the relation: 



CCF 



^Z,('-^)VCCFl 



Si) 



^JZlCCF(i) 



scale : 



(A.2) 



where CCF(i) is the Cross-Correlation Function measured for 
the velocity / and CCF„oise{i) is the quadratic sum of the pho- 
ton noise and the detector noise integrated inside the CCF mask 
holes for the velocity /. The correction factor N scale corresponds 
to the scale of the velocity step in detector pixel unit. The RV 
uncertainty 6Vrms is then equal to: 



6Vr 



1 



Qccf VS CCF{i) 



(A.3) 



This RV uncertainty, directly computed on the CCF, is ro- 
bust and can be applied to any kind of spectra down to SNR of 
~20. Figure IATTI shows the RV uncertainty scaled to an equiva- 
lent SNR of 100 and computed for a sub-sample of our planet 
search program with SOPHIE in HR mode. The RV uncertainty 
is plotted versus the index B-V and for two classes of v sin / (<3 



«1 
I' 



S/N- 100 



0,8 
B-V 



Fig. A.l. RV uncertainty versus the index B-V and for two 
classes of vsin/, <3 kms ' (plain dark circles) and in between 
6 and 10 kms ' (open blue circles), for SOPHIE observations 
in High Resolution mode with SNR of 100. 

km s ' and in between 6 and 10 km s~'). It illustrates the depen- 
dance with the effective temperature and the projected rotational 
velocity of the stars. 

Appendix B: Determination of vsin / and [Fe/IH] from 
tlie SOPHIE Cross-Correlation Function 

We used the property that the CCF can be considered as an aver- 
age line of the target spectrum. Line profiles are directly related 
to the stellar atmospheric parameters (e.g. abundances, thermal 
expansion), macro and micro-turbulence, rotation rate and in- 
strumental profile. Then, stellar properties affect spectral lines 
included in the correlation mask and are reflected in the CCF. 
We can derive stellar properties in analyzing the FWHM and the 
contrast of the fitted Gaussian on the SOPHIE CCF. The v sin / 
of a star might be derived as described by Benz & Mayor (1981) 
for the CORAVEL spectrograph. They refined their calibration 
afterward (Benz & Mayor 1984). Queloz et al. (1998a) applied it 
to the ELODIE spectrograph. The stellar [Fe/H] calibration was 
first done by Mayor ( 1 980) for the CORAVEL CCF, improved by 
Pont (1997). The both calibrations is done here for the SOPHIE 
spectrograph following the Santos et al. (2002) methodology for 
the CORALIE spectrograph. The spectral line profile is depen- 
dent of the instrumental one, so these calibrations had to be de- 
rived for diff'erent instruments. 



B.1. Calibration of tine projected rotational velocity v sin I 

We calibrated the relation between the vsin/ and the width cr 
of the CCF for both observing modes. High Efficiency (HE) 
and High Resolution (HR), of SOPHIE spectrograph using the 
cross-correlation masks G2 and K5 which have almost the same 
CCF width (negligible at the photon noise level). We note that 
the metallicity does not affect too much the width of the lines 
in comparison the the expected precision of our calibration. For 
solar-type stars, these two variables can be related by: 



V sini= A 



yjcr^ - crl (B.l) 

where cr represents the measured Gaussian width of the CCF 
(FWHM/2 V21n2), ctq the value of the expected cr for a "non- 
rotator" of a given spectral type, and A is a constant relating the 
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Fig. B.l. Gaussian width of the SOPHIE CCF versus B-V index 
for 240 dwarf stars observed in HR mode. The lower envelope 
determines the locus of "non-rotator" stars. 
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Fig. B.2. Known v sin / versus measured cr in the HR mode and 
the adjusted slope. 

"excess" width of lines to the actual projected rotational velocity 
V sin i. 

We first determine the quantity ctq by adjusting the lower 
envelope of the distribution of points in the diagram, shown 
in Figure IB. II for HR mode, of cr versus B-V which is quite 
well related with stellar temperature for dwarfs with B-V 
index in the range 0.4 - 1.3. We neglected here the effect of the 
metallicity. The adjusted lower envelopes for HR and HE mode 
are well described by; 

o-Q_HR = 9.90-22.56(B-y)H-22.37(B-y)^-6.95(B-y)\B.2) 
o-o_HE = 10.52-22.56(B-y)-H22.37(B-y)^-6.95(B-y)^(B.3) 

The other variable in Eq. lB.ll we need to determine is A, i.e., 
the constant relating the v sin / to the excess width of the CCF. 
For that we observed with SOPHIE 14 stars with known v sin / 
from the catalog of rotational velocities of F and G stars (Reiners 
& Schmitt 2003) and we find A = 1.73 and A = 1.64 for respec- 
tively the HR and HE mode. Figure lB72l shows the relation ob- 
tained between the known v sin / and the measured cr for the HR 
mode. We note that this relation is valid up to vsin;=20kms ' 
(Bouvier, private communication). The estimated uncertainty on 
vsin/ is 1 kms"'. 

We find only three values of v sin / from known planet host- 
ing stars: HD 189733 and HD80606 from Butler et al. (2006) and 



HD17156 from Fischer et al. (2007) with spectroscopic determi- 
nation that may have its own uncertainty. Their estimations for 
these stars are about 1.5 kms"' above our estimate. On the other 
hand, our calibration find a value of 2.6 kms"' for the Sun v sin /. 
These values are satisfactory for the precision required. 

The technique described above represents a simple way of 
obtaining an estimation of projected rotational velocities for 
dwarfs, simply as a by-product of the precise radial-velocity 
measurements. In particular, this calibration permits us to ob- 
tain quite easily values for the vsin / for all the stars in the 
SOPHIE exoplanet search program. 

B.2. Calibration of tiie metallicity 

The different SOPHIE masks are mainly built with neutral weak 
lines. An important quantity of ferric ions increases the absorp- 
tion coeflicient for the lines of these elements. The lines then 
become deeper and the contrast of the CCF larger. But, the line 
contrast depends also on the spectral type, measured by the B-V 
for a same luminosity class. The higher is the B-V, the lower 
is the effective temperature and the spectrum has more lines of 
metallic atoms. In addition, stellar rotation and Doppler shift, 
related to the relative speed of each atoms, modify the shape 
of the CCF widening and spreading the line profile but, it does 
not change the line surface. To take into account the stellar ro- 
tation and the macroturbulence at the stellar surface, we com- 
puted the total area W of the Gaussian in spite of CCF con- 
trast, W - V27r(contrast/100)cr where cr represents the mea- 
sured Gaussian width of the CCF and the contrast is expressed 
in percentage, and studied its dependance as a function of the 
B-V and metallicity. It is expected that the CCF surface grows 
with metallicity and with the B-V. 

A unique independent measurements of the metallicity has 
to be done for the calibration. We used a spectroscopic analy- 
sis based on a detailed study of each line (Santos et al. 2004). 
We determined accurately the [Fe/H] of 32 targets covering the 
B-V and CCF surface domains. The mean uncertainty from the 
spectroscopic determination is about 0. 1 dex. 

SOPHIE has two modes of observations (Bouchy et al. 
2009b), the High Resolution (HR) and the High Efficiency (HE) 
mode. Two fibers are used for each mode: one for the star and 
the other for the sky spectrum or for simultaneous calibration 
lamp exposure. The CCF surface is independent of the observa- 
tion mode, provided that the SNR is sufficient (i.e. SNR > 20 or 
SNR > 100 for thorium simultaneous mode which produces an 
important background of scattered light all over the CCD). On 
the other hand, the CCF contrast, and so the CCF surface, de- 
pends on the calibration mask used. We studied SOPHIE masks 
of spectral type G2 and K5 because the calibration is limited in 
B-V, excluding the coldest and the warmest stars. For the K5 
mask, the fit, shown in Fig. IB. 31 gives: 

[Fe/H] = 0.2615 + 3.9553 log W - 2.2606(B - V) (B.4) 
Samely, for the G2 mask, we obtained: 

[Fe/H] = -0.9440 H- 3.8807 log ly- IA992(B-V) (B.5) 

These calibrations are valid for 0.43 <(B-V)< 0.98 and -0.47 
< [Fe/H] < 0.44. 

We compare the [Fe/H] values obtained from the CCF sur- 
face and from the spectroscopic analysis in Fig. IB. 41 The disper- 
sion of the residuals around the zero values is 0.09 dex which is 
satisfactory given the spectroscopic uncertainty. We do not ob- 
serve any trend with the stellar temperature or B - y as shown 
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in Fig. IB .41 Our calibration is in agreement with others spec- 
troscopic determination of the metallicity: Sousa et al. (2006), 
Sousa et al. (2008), Gonzalez et al. (2007) and Fischer & Valenti 
(2005) as shown in Fig.lO] 

The tool included in the SOPHIE data reduction software al- 
lows an immediate estimation of the [Fe/H] of the star with a 
good approximation (± 0.09 dex). Our aim is to have a quick es- 
timate of the target metalUcity to focus our planet-search survey 
around over-metallic stars. One goal of the program is to find 
new hot-Jupiters known to be more prevalent around metal-rich 
stars (e.g., Santos et al. 2004; Santos et al. 2005). 



Appendix C: Determination of log/?^,^ from the 
SOPHIE spectra 

The SOPHIE spectra, with a bandpass from about 3900 A to 
6800 A, include both Calcium II H and K resonant lines cen- 
tered at 3968.49 A and 3933.68 A respectively. These lines are 
widely used indicators of stellar magnetic activity. In Fig. IC. H are 
plotted two high-SNR SOPHIE specti-a of the Ca II lines central 
region for the chromospherically active star HD 131156 and the 
inactive HD 187013. An active star for which spots and plages 
are present on the photosphere is characterized by an emission 
in the center of the absorption line (e.g. Boisse et al. 2009). 

The stellar activity is parametrized by an index created by 
Wilson (1968) defined by the ratio between the emitted flux in 
the center of the lines and the continuum flux. The S-value is 
defined by the measure of the quotient of the flux in two trian- 
gular bandpasses centered on the H and K emission cores and 
two continuum regions on either side. We define our Ca II H&K 




0.4 0.6 0.8 1 



B-V 

Fig.B.3. Surface of the CCF, W, as a function of the B - V 
for the K5 correlation mask. The dot size is proportional to the 
spectroscopic metallicity. The solid lines draw the relationship 
between the parameters for [Fe/H]=cst. The fit with an exponen- 
tial is good in the domain of metallicity but it is obviously wrong 
forB - y>1.0. 
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Fig.B.4. Top: [Fe/H] computed from the calibration as a func- 
tion of the spectroscopic estimation. The solid line illustrates a 
1 : 1 correlation. Bottom: Residuals from the calibration as a func- 
tion of the B-Y. 

index following the Mt. Wilson Smw index (Baliunas et al. 1995; 
Boisse et al. 2009): 

H -t- ^ 

Index — (C.l) 

B + y 

where H and K are the flux measured in the triangular 1.09 A 
FWHM window centered on each line of the Ca II doublet, and 
B and V estimate the continuum on both sides with the flux mea- 
sured in 20 A wide windows respectively centered on 3900 A 
and 4000 A. An automatic cut of the cosmic rays are done. Each 
line is located in two consecutive orders of the SOPHIE spectra. 
For the K line, we average the flux measured in two orders. As 
one of the H line is located on the edge of the an order where less 
flux is collected which introduce noise, only the flux measured 
in the H line of another order is kept. 

When studying the Ca II H&K lines in SOPHIE spectra, one 
must account for the fact that they are located in a spectral region 
that suffers from contamination from background scattered light. 
In order to use the signature of flux emitted in the Ca II lines, we 
estimate a background light level and subtract it from the stellar 
spectra. For each mode, two fibers are used on SOPHIE, one of 
which is pointed on the target. The other fiber receives the sky 
spectrum or is fed with a Thorium- Argon lamp. The background 
light level of a stellar order is estimated from the same order of 
the other fiber by fitted a polynomial function on local minima. 
We put limits on the minimal SNR in the first order (/l~3955A) 
where both Ca II H&K lines are found to have a good estimation 
of the activity index. The boundaries are fixed where dependency 
between SNR and log is no longer observed. If the second 
aperture received the sky spectrum, the background light level 
is low and the minimal SNR is equal to 10. When the Thorium- 
Argon lamp feed the second fiber, the level of background light 
on the CCD is about 0.1% of the stellar continuum and the min- 
imal SNR to have a reliable value of log R'^^y^ is fixed at 30. 

We measure a chromospheric index S sophie that is calibrated 
to the Mount Wilson index S'mw by observing calibration stars 
in HE and HR SOPHIE mode. The Mount Wilson values from 
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Fig. B.5. Top: [Fe/H] from others estimations from spectro- 
scopic analyses: Sousa et al. (2006) (black squares), Sousa et al. 
(2008) (red open triangles), Gonzalez et al. (2007) (blue circle) 
and Fischer & Valenti (2005) (filled green triangles) as function 
of those from SOPHIE calibration. The solid line illustrates a 1 : 1 
correlation. Bottom: Residuals from the top panel as a function 
of [Fe/H] from SOPHIE calibration. The rms of the residuals is 
0.09 dex. 




Fig. C.2. Mean 5sophie values as function of the Mount Wilson 
values from Duncan et al. (1991). The solid line is the best-linear 
fit to the data. Left: High Resolution mode, (rms - 0.012) Right: 
High Efliciency mode, (rms = 0.007) 

Duncan et al. (1991) and our measurements 5 sophie are listed in 
Table lC.2l available in the electronic form of the paper Table lC.2l 
contains in its cols. 1-6, the name of the stars, the S'mw value and 
its error, the S sophie and its error and the SOPHIE observational 
mode, respectively.. The mean 5 sophie values compared with 
5mw are shown in Fig. IC.2l for the set of calibrating stars for 
both mode of SOPHIE. In HR mode, the best linear fit to the 
data gives: 

5 SOPHIE = (0.71 + 0.03)5 Mw + (0.039 ± 0.006) (C.2) 
In HE mode, the relation is: 

^SOPHIE = (0.72 + 0.02)5mw + (0.021 ± 0.006) (C.3) 

With this index, we can calculate the litterature's index cor- 
rected from photospheric emission, R'^^^ (Noyes et al. 1984). 
This measure is available only for 0.43 < (B - V) < 1.2. 

We estimate the uncertainty of our index (±0. 1 dex) with the 
dispersion on our values for a same star over one year This dis- 
persion may be due to intrinsic variations of the stellar activity. 
But, we consider that our main error comes from the subtraction 
of the background light. 




3920 



3940 3960 
Wavelength [A] 



3980 



Fig. C.l. Two SOPHIE spectra of the region of the Call H and 
K lines between 3900 and 4000 A. Bottom: Active star. Top: Non 
active star 
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Table C.l. Radial velocities of HD 109 246 measured with SOPHIE. 



BJD RV ±lo- 

-2400000 (kms-') (kms-') 

54126.66704 -19.4950 0.0064 

54127.61475 -19.4953 0.0061 

54506.61615 -19.4215 0.0063 

54525.49634 -19.4682 0.0063 

54525.63339 -19.4727 0.0057 

54536.58945 -19.4923 0.0059 

54545.51957 -19.4961 0.0057 

54554.47615 -19.4745 0.0057 

54555.49364 -19.4929 0.0054 

54584.45332 -19.4474 0.0055 

54597.46134 -19.4810 0.0059 

54642.41898 -19.4347 0.0057 

54646.36523 -19.4415 0.0057 

54646.37998 -19.4597 0.0056 

54662.35377 -19.4669 0.0058 

54664.35273 -19.4551 0.0066 

54666.38121 -19.4605 0.0059 

54697.32306 -19.4539 0.0058 

54717.29649 -19.4247 0.0057 

54722.28695 -19.4175 0.0067 

54724.28200 -19.4404 0.0059 

54822.66825 -19.4878 0.0059 

54834.69504 -19.4464 0.0059 

54834.72295 -19.4436 0.0059 

54835.70783 -19.4411 0.0059 

54852.62557 -19.4238 0.0064 

54853.69041 -19.4287 0.0060 

54853.71034 -19.4312 0.0061 

54879.60628 -19.5014 0.0059 

54880.58603 -19.4929 0.0060 

54882.61111 -19.5078 0.0059 

54883.69930 -19.5128 0.0058 

54884.55744 -19.5085 0.0059 

54886.47681 -19.4934 0.0068 

54887.52148 -19.5006 0.0064 

54888.58770 -19.4891 0.0059 

54889.61323 -19.4865 0.0059 

54890.61440 -19.4955 0.0058 

54893.48340 -19.4834 0.0059 

54902.50858 -19.4430 0.0059 

54903.54493 -19.4270 0.0059 

54904.47549 -19.4412 0.0059 

54905.48583 -19.4326 0.0058 

54906.52068 -19.4311 0.0059 

54911.52222 -19.4232 0.0063 

54912.49518 -19.4202 0.0063 

54924.57038 -19.4315 0.0065 

54926.45707 -19.4679 0.0056 

54934.48378 -19.4848 0.0059 

54941.55035 -19.4873 0.0058 

54946.35550 -19.4947 0.0059 

54968.38892 -19.4680 0.0057 

54971.38402 -19.4588 0.0057 

55004.36914 -19.4659 0.0060 

55029.35623 -19.5058 0.0056 

55238.54345 -19.4777 0.0076 

55240.69127 -19.4646 0.0059 

55242.61572 -19.4487 0.0059 
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Table C.2. Stars used to calibrate the Ssophie in S Mount wuson- 



Stars 


S MW 


Mw) 


S SOPHIE 


Sophie) 


Mode 


HD 97334 


0.328 


0.012 


0.276 


0.013 


HR 


HD 101501 


0.310 


0.021 


0.264 


0.021 


HR 


HD 114378 


0.238 


0.007 


0.185 


0.007 


HR 


HD 114710 


0.200 


0.023 


0.198 


0.023 


HR 


HD 115404 


0.546 


0.005 


0.420 


0.005 


HR 


HD 126053 


0.168 


0.004 


0.179 


0.004 


HR 


HD 131156 


0.449 


0.022 


0.362 


0.016 


HR 


HD 141004 


0.159 


0.006 


0.153 


0.003 


HR 


HD 143761 


0.148 


0.001 


0.141 


0.027 


HR 


HD 149661 


0.356 


0.041 


0.320 


0.004 


HR 


HD 152391 


0.392 


0.029 


0.349 


0.006 


HR 


HD 154417 


0.271 


0.010 


0.266 


0.010 


HR 


HD 182101 


0.212 


0.005 


0.192 


0.001 


HR 


HD 187013 


0.150 


0.002 


0.138 


0.002 


HR 


HD 190007 


0.775 


0.067 


0.592 


0.070 


HR 


HD 194012 


0.195 


0.004 


0.185 


0.004 


HR 


HD 201091 


0.617 


0.065 


0.465 


0.070 


HR 


HD 101501 


0.310 


0.021 


0.254 


0.021 


HE 


HD 114710 


0.200 


0.023 


0.158 


0.101 


HE 


HD 115404 


0.546 


0.005 


0.409 


0.005 


HE 


HD 131156 


0.449 


0.022 


0.363 


0.003 


HE 


HD 141004 


0.159 


0.006 


0.133 


0.002 


HE 


HD 187013 


0.150 


0.002 


0.125 


0.002 


HE 


HD 194012 


0.195 


0.004 


0.173 


0.004 


HE 



